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We consider a trapped atomic Bose-Einstein condensate penetrated by a repulsive Gaussian 
potential and theoretically investigate the dynamics induced by oscillating the Gaussian potential. 
Our study is based on the numerical calculation of the two-dimensional Gross-Pitaevskii equation. 
Our calculation reveals the dependence of the characteristic behavior of the condensate on the 
amplitude and frequency of the oscillating potential. These dynamics are deeply related to the 
nucleation and dynamics of quantized vortices and solitons. When the potential oscillates with a 
large amplitude, it nucleates many vortex pairs that move away from the potential. When the 
amplitude of the oscillation is small, it nucleates solitons through annihilation of vortex pairs. We 
discuss three issues concerning the nucleation of vortices. The first is the phase diagram for the 
nucleation of vortices and solitons near the oscillating potential. The second is the mechanism and 
critical velocity of the nucleation. The critical velocity of the nucleation is an important issue in 
quantum fluids, and we propose a new expression for the velocity containing both the coherence 
length and the size of the potential. The third is the divergence of the nucleation time, which is the 
time it takes for the potential to nucleate vortices, near the critical parameters for vortex nucleation. 
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I. INTRODUCTION 

The Gross-Pitaevskii (GP) equation has a long his- 
tory in the study of quantized vortices and solitons, and 
is important for both condensed matter and nonlinear 
physics. In superfluid ^He and ^He, vortices and solitons 
appear and numerous studies have been made of them 
[ll, [ij]. However, superfluid ^He, which exhibits strong 
interactions between particles, cannot be quantitatively 
described by the GP equation since it can be applied only 
to dilute Bose systems. 

A dilute atomic Bose-Einstein condensate (BEC) has 
been experimentally realized [3.-^] and is quantitatively 
described by the GP equation. Therefore, the dynam- 
ics of vortices and solitons expected by the GP equation 
can be experimentally observed. Thus, the experimental 
realization of dilute atomic BECs sets the stage for the 
theoretical study of the GP equation. 

Vortices in atomic BECs have been experimentally and 
theoretically studied. The nucleation and dynamics of 
vortices are mainly treated in two systems. One is a sys- 
tem where the trapping potential rotates and in the other 
a localized potential moves in the atomic BEC. In the for- 
mer case, the vortices are nucleated from the surface of 
the condensate, and the vortex lattice is formed through 
the characteristic nonlinear dynamics [6, 7]. In the latter 
case, vortexpairs are nucleated by the uniform moving 
potential [1,[9|], giving rise to interesting dynamics. Neely 
et al. [lo( experimentally and numerically observed that 
the vortices nucleated by the potential migrate with long 
periodicity in the condensate. The Karman vortex street 
in atomic BECs has been studied by numerical calcula- 
tions based on the GP equation pj| . 

Solitons have also been investigated in atomic BECs 
[l2|-0^] . Solitons of the GP equation in a one-dimensional 



system are stable, while those in two or three-dimensional 
systems are unstable against long wavelen gth transverse 
oscillation, referred to as sneak instabilitv|l2l, [l5| . The 
instability of solitons leads to the nucleation of vor- 
tices, which has been experimentally and numerically 
confirmed [H, [l3]- Thus, solitons are profoundly related 
to vortices. 

The solitary wave solution of the GP equation has been 
investigated in detail by Jones and Roberts [16]. They 
calculated the energy Es and the momentum Pg of the 
solution, and found that there are two branches which 
are drawn from a common point in the Eg-Pg plane, the 
upper and lower branches, which denote the rarefaction 
pulse and the vortex pair, respectively pj|, [l7|. We ex- 
pect that a vortex pair can dynamically change into a 
rarefaction pulse and vice versa. We have numerically 
observed that annihilation of a vortex pair creates a rar- 
efaction pulse and the collapse of a pulse leads to the 
nucleation of a vortex pair, which is caused by an oscil- 
lating potential [19]. 

A currently important theme in quantum hydrody- 
namics is quantum turbulence (QT), which was first in- 
vestigated in superfluid helium [20] . There are a number 
of methods to create QT in this system. One is to employ 
thermal counterflow. Another is to oscillate objects such 
as grids, micro spheres, wires, and forks in superfluid 
helium. Rece ntly , QT in an atomic BEC has been vigor- 
ously studied [2ll-[23|. There are also several methods of 
creating QT in this system, namely phase printing plj , 
precessing a trapping potential [22], and oscillating the 
potential [23]. To date, experimental creation and obser- 
vation of QT in atomic BECs has been achieved only by 
oscillating a trapping potential [23]. 

We apply a method for dealing with oscillating objects 
in superfluid helium to atomic BECs. However, there is a 
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definite difference between atomic BECs and helium sys- 
tems. In helium, moving objects do not nucleate vortices 
but simply amplify remanent vortices which already exist 
in the system, creating QT [Ts*]. It is difficult to theo- 
retically treat the nucleation of vortices in helium due to 
the strong correlation between the atoms and it is not 
easy to control the event experimentally [1]. In contrast, 
moving objects intrinsically nucleate vortices in experi- 
ments with atomic BECs and optical techniques enable 
us to observe them [lo[. The nucleation of vortices and 
their dynamics are quantitatively described by the GP 
equation. In this system, a blue detuned laser is used 
as the object, which is represented by a repulsive Gaus- 
sian potential in our numerical calculation. We are in- 
terested in the nucleation and dynamics of vortices and 
solitons induced by oscillations of the potential, which 
should depend on the amplitude and frequency of the 
oscillating potential. Previously, we have reported that 
the oscillating potential causes synergy dynamics of vor- 
tices and solitons peculiar to the oscillation, though the 
dependence on the parameters has not been investigated 
flQj . In this paper, through the investigation of the dy- 
namics induced by the oscillating potential, we study the 
dependence of the dynamics on the parameters and also 
the mechanism for the nucleation of vortices, the critical 
velocity of the nucleation, and the divergence of the nu- 
cleation time, which have not been satisfactorily treated 
for atomic BECs. 

Previous research has been made on the response of 
atomic BECs induced by an oscillating potential [25l - [27j . 
However, the previous studies focused on dissipation and 
drag forces and barely considered the instability and non- 
linear dynamics induced by an oscillating potential, in- 
cluding dynamical instability. Landau instability, nucle- 
ation and dynamics of vortices and solitons, and quantum 
turbulence. Therefore, our study is very different from 
previous studies. 

The nucleation of a vortex pair by an oscillating poten- 
tial can be related to dynamical critical phenomena. In 
fact, nucleation by a linear moving potential has been in- 
vestigated in terms of the phenomenon, and a power law 
between the vortex emission frequency and the velocity 
of the object near the critical velocity has been reported 
through a numerical simulation of the GP equation [28|. 
Thus, we address the nucleation by an oscillating poten- 
tial in terms of dynamical critical phenomena. Moreover, 
turbulence is deeply related to pattern formation [29], 
so that the dynamics induced by the oscillating poten- 
tial is clearly connected with the process. These phe- 
nomena have been historically investigated by the com- 
plex Ginzburg-Landau equation, whose typical example 
is the GP equation. Thus, it is possible to study non- 
linear and nonequilibrium phenomena, namely pattern 
formation and dynamical critical phenomena, in atomic 
BECs. Therefore, our study can pioneer a new direction 
in atomic BECs. 

The nucleation of vortices in quantum ffuids has been 
considered to be an important issue, including the mech- 



anism and critical velocity of the nucleation. Previously, 
nucleation has been investigated in superffuid ^He, where 
the core size of the vortices is so small that their visu- 
alization is hard. The theoretical study of nucleation is 
also difficult since the interaction between the "^He atoms 
is very strong. However, in atomic BECs, this situation 
changes. Optical techniques make visualization possible 
and the dynamics of vortices can be followed in experi- 
ments. The nucleation and dynamics of vortices are eas- 
ily treated because atomic BECs are well described by 
the GP equation. In this paper, we focus on the mecha- 
nism and critical velocity of nucleation. The mechanism 
is discussed in terms of the divergence of the quantum 
pressure. Our numerical calculations reveal that the co- 
herence length and size of the potential are very relevant 
to the critical velocity of the nucleation, and we propose 
an expression of the velocity including both lengths. 

We consider an atomic BEG confined by a harmonic 
trap to investigate the dynamics of the condensate in- 
duced by a repulsive oscillating Gaussian potential. Our 
study assumes that the condensate is pancake shaped 
and the temperature of the system is nearly zero, which 
means that the condensate is well described by the two- 
dimensional GP equation. 

This paper is organized as follows. In Sec. II, we de- 
scribe the formulation. Section HI treats the breakdown 
of Kelvin's theorem on circulation, which leads to the 
nucleation of vortices. In Sec. IV, the synergy dynamics 
of vortices and solitons are described. The nucleation of 
solitons near the oscillating potential is discussed in Sec. 
V. Section VI describes the nucleation of multiple vortex 
pairs. A phase diagram summarizing Sees. IV- VI is pre- 
sented in Sec. VII. The critical velocity for the nucleation 
is treated in Sec. VIII. In Sec. IX, we show the results 
on divergence of the nucleation time. Section X discusses 
the heating of the condensate. The dependence of the as- 
pect ratio of the condensate on the dynamics, the noise 
of the system, and the appropriate parameters for the 
growth of QT are discussed in Sec. XI. We summarize 
our study in Sec. XIL 



II. FORMULATION 

We consider an atomic BEG near zero temperature, 
in which the interaction between particles is very weak 
because of the low density and is replaced by an effec- 
tive interaction proportional to the delta function and 
the s-wave scattering length. Almost all particles are 
condensed near zero temperature so that this system is 
well described by a macroscopic wavefunction obeying 
the GP equation. In this work, we assume that the con- 
densate is strongly confined along the z direction, that 
is, a quasi two-dimensional system. Therefore it is well 
described by the two-dimensional GP equation: 

'^? = -T^^'^ + ^^ + 5l^lV, (1) 
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FIG. 1: (Color online) Initial density profile is shown in (a), 
where the x and y dimensions of the images are 145 /xm and 
34.0 /im, respectively, (b) shows the enlarged figure of the 
square box in (a), with x and y dimensions of 18.0 /im and 10.0 
/im, respectively. The Thomas-Fermi radii Rx and Ry are 
69.3 /im and 13.9 /xm, respectively. The low density region in 
the center of the condensate is formed by a repulsive Gaussian 
potential. The diameter of the circle in (b) is 5 /xm. 



where m is the particle mass, V is the potential, and g 
is an interaction parameter for the two-dimensional case, 
where the wavefunction along the z direction is assumed 
to be Gaussian [7|]. The wavefunction i/j is normalized 
by the total particle number N. We suppose that the 
condensate is confined by a harmonic potential Vh and 
penetrated by a Gaussian potential Vq , so that V = Vh-\- 
Vg where Vh = \m{ijj'^x^ -\- oUyy'^) and 



Vg = Voexp[-{{x - xo{t))' + y'}/£]. 



(2) 



Here xo{t) is the x-coordinate of the center of the Gaus- 
sian potential and d is its radius. We set the Gaussian 
potential to oscillate as xo{t) = esin(cc;^). We define the 
velocity of the oscillation as v = ecj. 

We set g = 4.05 x IQ-"^^ J/m^ m = 1.42 x IQ-^^kg, 
iV = 8.0 X 10"^, cj^ = 27r X 4 /s, ujy = 27r x 20 /s, 
d = 1.30 /im, and Vq = lO^no as parameters. Here no 
is the density near the center of the condensate. We use 
a dimensionless form of Eq. (1) to perform a numeri- 
cal calculation with the Crank-Nicholson method. Space 
and time are normalized by h/ ^2mgnQ and h/gno^ and 
the space in the x and y directions is discretized into 
2048x640 bins. 

Figure 1 shows the initial density profile in which the 
low density region of the center is formed by the repul- 
sive potential. This is the ground state before the po- 
tential starts to move, obtained by the imaginary time 
step method. The size of the potential L is defined as 
the rounded value of the diameter of the boundary of 
the low density region, which is lower than half of the 
bulk density. In our case, L = 5 /im. Figure 1(b) shows 
that the potential of size L is surrounded with a bound- 
ary layer, which is a low density region with a width of 
the order ^. This boundary layer plays a key role for 
any dynamics of vortices near the potential, as discussed 
later. 

The accuracy of the numerical calculation is confirmed 



by the two following methods. One is the mirror symme- 
try about the x axis, which must be maintained because 
the initial state has mirror symmetry and the potential 
oscillates in the x-direction to maintain the symmetry. 
The other is to check that the final state of the conden- 
sate, which is obtained by the numerical calculation with 
forward time evolution, returns properly to the initial 
state in Fig. 1 when we perform the calculation with 
backward time evolution from the final state. This is 
confirmed by the total energy and the momentum in the 
X direction. 



III. BREAKDOWN OF KELVIN'S THEOREM 
ON CIRCULATION 

Vortices can be nucleated in our system because of the 
breakdown of Kelvin's theorem on circulation. This the- 
orem is established in classical fluid dynamics described 
by the Euler equation and ensures the conservation of 
vorticity, which leads to neither nucleation nor the dis- 
appearance of vortices. The breakdown is caused by the 
divergence of the quantum pressure of the GP equation 
as shown in the following. 

First, we derive the general expression for the Lagrange 
derivative of the circulation in the quantum hydrodynam- 
ics described by the GP equation. Using the Madelung 
transformation ijj = ^/nex.^p{i<p) where n and (j) are the 
density and the phase of the wavefunction, we obtain 
from Eq. (1) the Euler-like equation 



with 



m— v = - V(/i + -mv ) 



(1 = V -\- ng 



2m^/n 



(3) 



(4) 



where the superfluid velocity v is given by JiV(j)/m. The 
characteristic term of the equation is the quantum pres- 
sure term V^-^/n, which does not appear in the 
Euler equation. Using these equations and the circula- 
tion n{C{t)) = J^^^^ V -dl, we can calculate the Lagrange 
derivative of the circulation 



f^{C{t)) = lim ^[ / v{r,t^St)-dl- [ v(r,t)'dl\, 

5t^Q dt Jc{t+5t) Jc{t) 



D 
Dt' 

(5) 

where C{t) is a closed contour, which rides the superfluid 
velocity fleld to change the conflguration, deflned by r = 
p{s^t) and dp{s^t)/dt = v{p{s^t)^t). Here 5 is a time- 
independent parameter denoting a point on and its 
change through the domain sq < s < Si with arbitrary 
constants sq and si represents the whole closed contour 
C{t) [30|. We rewrite Eq. (5) with the parameter s 
to transform the integration range of the two terms to 
the same range of the domain sq < s < si, which make 
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possible the subtraction in Eq. (5). As a result, we obtain 
the following expression: 

^i^(C) = -— I V(jl)-dl^ [ {(Vxv)xv}-dl^ [ V'dv. 
Dt mJc Jc Jc 

(6) 

This is the general expression for the Lagrange deriva- 
tive of the circulation in the quantum hydrodynamics 
described by the GP equation. 

Second we consider the condition of the nucleation of 
vortices with Eq. (6), which is derived by the breakdown 
of Kelvin's theorem on circulation. Now, we treat the 
system without vortices to consider the condition of the 
nucleation of vortices. In the system, the second and 
third terms in the right hand of Eq. (6) vanish since the 
velocity v is potential flow and a single-valued function 
of the position. Unless a low density region causing the 
divergence of the quantum pressure is formed at a point 
along the contour C, jl in Eq. (6) does not diverge on 
the C, which leads to Dn/Dt = 0. In this case, vortices 
cannot be nucleated because of Kelvin's theorem on cir- 
culation. However, Dn/Dt does not generally vanish if a 
low density region causing the divergence of the quantum 
pressure is formed on C. This allows nucleation of the 
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FIG. 3: (Color online) Reconnect ion of vortex pairs near the 
oscillating potential: The density (left) and phase (right) pro- 
files at (a) t = 31.4 ms, (b) t — 38.3 ms , and (c) t — 43.2 ms 
are shown, where the x and y dimensions of each image are 
27.0 /im and 14.0/im. The symbols — and + denote a vortex 
with clockwise or counterclockwise circulation, respectively. 
The black arrows indicate the direction of motion of the po- 
tential. The value of the phase varies from — tt (white) to tt 
(black). The density and phase profiles before the collision 
between the potential and the vortex pair are shown in (a). 
Thereafter, another ghost vortex pair nucleates in (b), exiting 
the potential through the collision, which causes reconnection 
of the vortices. As a result, two pairs appear in (c). 



vortices, but the divergence is only the necessary condi- 
tion for the nucleation of the vortices (3l|. 

This situation is directly described in our numerical 
simulation. The low density region in the condensate ap- 
pears through two patterns. One is that the oscillating 
potential forms the low density region because of repul- 
sion, which causes nucleation of vortices. The other is 
the inside of the solitons, which allows the transforma- 
tion between vortices and solitons. 



IV. SYNERGY DYNAMICS OF VORTICES 
AND SOLITONS 



FIG. 2: (Color online) Nucleation of a vortex pair by the oscil- 
lating potential: The density (left) and phase (right) profiles 
at (a) t — 4.71 ms, (b) i — 8.05 ms , and (c) t — 13.7 ms are 
shown, where the x and y dimensions of each image are 23.0 
/im and 14.0 /xm. The symbols — and + denote a vortex with 
clockwise or counterclockwise circulation, respectively. The 
black arrows indicate the direction of motion of the potential. 
The value of the phase varies from — tt (white) to tt (black). 
A ghost vortex pair nucleates inside the potential (a), exits 
it (b), and finally fully leaves the potential (c). The parame- 
ters of the oscillation are e = lO/xm and cj = 60/s, which are 
common through Figs. 2-8. 



An oscillating potential creates vortex pairs, causes re- 
connection of pairs characterized by the oscillation, and 
causes the new pairs to leave for the surface of the con- 
densate. Consequently, the surface becomes filled with 
vortices having positive and negative circulation, which 
leads to nucleation of solitons and the migration of vor- 
tices. We call this sequence synergy dynamics of vortices 
and solitons, which often occurs in cases where the ampli- 
tude of the oscillation is larger than L. We demonstrate 
these dynamics through a simulation with e = 10/im and 
uj = 60/s. The overall dynamics have been briefly re- 
ported in [Toj, but this section describes each process in 
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FIG. 4: (Color online) Separation of vortex pairs near the sur- 
face of the condensate: The density profiles at (a) t = 58.9 ms 
and (b) t = 88.4 ms are shown, where the x and y dimensions 
of each image are 145 /xm and 34.0 /im, respectively. The 
symbols — and + denote a vortex with clockwise or counter- 
clockwise circulation, respectively. The black arrows indicate 
the direction of motion of the vortices. The vortex pairs ap- 
proach the surface of the condensate in (a) and separate and 
leave for the bow of the condensate in (b). 



the excitation of the surface waves, leading to the cre- 
ation of usual vortices [gI, Thus, the periphery of the 
condensate provides seeds of topological defects. In our 
system, the oscillating potential provides seeds within it- 
self. The potential starts to move, inducing a velocity 
field like back-flow, emitting phonons, and a ghost vor- 
tex pair is nucleated inside the potential as shown in Fig. 
2(a). The ghost pair tends to move away from the poten- 
tial in Fig. 2(b), and a usual vortex pair appears in the 
condensate in Fig. 2(c). Thus, the ghost vortices work 
as seeds of usual vortices. 

Reconnection of vortex pairs occurs near the oscillat- 
ing potential. The new vortex pair has an impulse in 
the same direction as that of the potential. Then, the 
potential changes the direction of the velocity. Thus, the 
potential will collide with the pair in Fig. 3(a). Then, 
a new ghost vortex pair is nucleated inside the potential 
whose impulse is opposite to that of the usual vortex pair. 



more detail. 

Ghost vortices, namely quantized vortices in a low den- 
sity region, are important for the nucleation of the usual 
vortices in the bulk density region since nucleation re- 
quires seeds of vortices. In rotating BECs, ghost vortices 
are nucleated outside the condensate, entering it through 
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FIG. 5: (Color online) Nucleation of solitons: The density 
and phase profiles at (a) t = 110 ms, (b) t = 128 ms, and 
(cl)-(c2) t = 138 ms are shown, where the x and y dimen- 
sions of each image are 45.2 /im and 27.6 /xm, respectively. 
(c2) is the phase profile corresponding to the density profile 
(cl). The symbols — and + denote a vortex with clockwise or 
counterclockwise circulation, respectively. The black arrows 
indicate the direction of motion of the vortices and solitons. 
The value of the phase varies from — tt (white) to tt (black). 
Some vortices sit near the surface in (a) and reconnection of 
the vortices occurs in (b), where the square boxes with dot- 
ted lines enclose the new vortex pairs. While the new pairs 
move toward the center of the condensate, the pairs annihi- 
late, which leads to nucleation of solitons in (cl) and (c2), 
shown by the square boxes with solid lines. 
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FIG. 6: (Color online) Structure of density and phase for nu- 
cleation of solitons: The density and phase profiles at (al)- 
(a3) t — 132 ms and (bl)-(b3) t — 135 ms are shown. The 
X and y dimensions of (al), (bl) are 27.0 /im and 22.0 //m, 
respectively, and for the rest are 13.0 /xm and 8.01 /xm, re- 
spectively. (a2) and (b2) show enlarged figures, depicted by 
square boxes with solid lines in (al) and (bl). (a3) and (b3) 
are the phase profiles corresponding to the density profiles 
(a2) and (b2). (al)-(a3) and (bl)-(b3) show the density and 
phase profiles before and after the annihilation of a vortex 
pair, respectively. The symbols — and + denote a vortex 
with clockwise or counterclockwise circulation, respectively. 
The black arrows indicate the direction of motion of the vor- 
tices and solitons. The value of the phase varies from — tt 
(white) to TT (black). (al)-(a3) show the density and phase 
profiles before the annihilation of vortices, and the profiles 
after the annihilation are shown in (bl)-(b3). 
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FIG. 7: (Color online) Collision and collapse of solitons: The 
density profiles at (a) t = 139 ms, (b) t = 147 ms , (c) t = 
157 ms, and (d) t = 167 ms are shown, where the x and y 
dimensions of each image are 45.2 /im and 27.6 /xm, respec- 
tively. The symbols — and + denote a vortex with clockwise 
or counterclockwise circulation, respectively. The black ar- 
rows indicate the direction of motion of the solitons. The 
square box with solid lines encloses the solitons and the box 
with dotted lines encloses the vortex pair, (a) and (b) show 
the state before and after the collision of the solitons, respec- 
tively. The decay of solitons to vortex pairs is shown in (c). 
(d) shows the density profile after the decay. 




FIG. 8: (Color online) Migration of vortices: The density pro- 
files at (a) t 98.2 ms, (b) t 236 ms , and (c) t 301 ms 
are shown, where the x and y dimensions of each image are 
145 /im and 34.0 /xm, respectively. The symbols — and + 
denote a vortex with clockwise or counterclockwise circula- 
tion, respectively. The black arrows indicate the direction of 
motion of the vortices. The vortices near the surface of the 
condensate in (a) move toward the bow, reconnecting there 
in (b) and returning to the center in (c). 



reconnecting with it as shown in Fig. 3(b). Thus, two 
new vortex pairs appear in the condensate in Fig. 3(c). 
Thereafter, the pairs move away from the potential, leav- 
ing for the surface of the condensate. This reconnection 
is characteristic of the oscillating potential because the 
potential repeatedly emits vortices of positive and nega- 
tive circulation in opposite directions, which is not seen 
for potentials of uniform motion. This leads to nucleation 
of solitons, as shown in the following. 

The vortex pairs separate as they approach the sur- 
face of the condensate in Fig. 4. This behavior is qual- 
itatively understood by applying the idea of an image 
vortex, which is often used in hydrodynamics. The vor- 
tices induce a circular velocity field in a uniform system, 
but the field is distorted in a nonuniform system. This ef- 
fect is strongly evident near the surface of the condensate 
where the density profile rapidly varies. As the vortices 
arrive at the surface, the normal component of the veloc- 
ity field is suppressed. This situation is approximately 
equal to the relation between a vortex and a solid wall, 
so that the dynamics of vortices near the surface in Fig. 
4 can be shown by the image vortex [32|. Note that this 
idea only gives a qualitative understanding since the sur- 
face is not exactly a solid wall. 

The vortices near the surface of the condensate have 
two fates. One is that a vortex pair transforms into soli- 
tons through the annihilation of the pairs. The other is 
that the vortices migrate in the condensate. We show 
these dynamics in the following. 

Transformation between vortices and solitons - Many 
vortices have accumulated near the surface of the conden- 
sate in Fig. 5(a) since the oscillating potential continues 



to make vortices with positive and negative circulation. 
Hence, the vortices near the surface can reconnect with 
each other as shown in Fig. 5(b), where we enclose the 
new vortex pairs with square dotted lines. These pairs 
have impulse toward the center of the condensate. As 
a pair approaches the center, the size of the pair dimin- 
ishes. Consequently, pair annihilation of vortices occurs, 
making the solitons shown in Figs. 5(cl) and (c2). 

A soliton of the GP equation with repulsive interaction 
has a locally decreased density profile and a rapidly vary- 
ing phase profile. The low density parts in Figs. 5(cl) 
and (c2) have these properties and hence we can identify 
them as solitons. This kind of nucleation of solitons is 
characteristic of an oscillating potential since it is caused 
by the potential emitting vortices in opposite directions. 

One may wonder why the annihilation of a vortex pair 
creates a soliton. Certainly, this would not occur if the 
system is uniform. However, our system is confined by a 
trapping potential, so that phenomena characteristic of a 
nonuniform system can occur, which is shown in Fig. 6. 
Figures 6 (al) and (bl) show the density profiles immedi- 
ately before and after the nucleation of solitons. We pay 
attention to the square boxes in Figs. 6(al) and (bl), 
which are enlarged in Figs. 6(a2) and (b2). The phase 
profiles corresponding to Figs. 6(a2) and (b2) are Figs. 
6(a3) and (b3), where the phase profiles are distorted 
because the vortices are in a narrow region between the 
high density region of the condensate and the surface. 
Fig. 6(a3) shows that the phase of the vortex pair in the 
front and rear is fiat while the phase in the right and 
left regions changes rapidly. In addition, the phase be- 
tween the pair also changes more rapidly than that in the 



right and left regions. Then, the velocity is suppressed in 
the front and rear, and enhanced in other regions, where 
the velocity between the pairs is larger than that in the 
right and left regions. Therefore, the vortices attract each 
other through the Magnus force, and annihilation of the 
pair occurs, leading to the distorted phase shown in Fig. 
6(b3) where the phase in the center of the low density 
region in Fig. 6(b2) rapidly changes by about 27r/3. 

After the nucleation of solitons, the solitons collide and 
collapse. As shown in Fig. 7(a), the solitons move toward 
the center of the condensate, so that collision occurs as 
shown in Fig. 7(a), (b), which show the density pro- 
file before and after the collision. The solitons do not 
change their configurations during the collisions, which 
is a property of solitons. Thereafter they move towards 
the surface of the condensate, decaying into new vortex 
pairs as shown in Fig. 7(c). At the surface, the pairs sep- 
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FIG. 9: (Color online) Nucleation of solitons near the oscil- 
lating potential: The density (left) and phase (right) profile 
at (a) t = 23.6 ms, (b) t = 27.5 ms , (c) t = 30.3 ms, and (d) 
t — 33.2 ms are shown, where the x and y dimensions of each 
image are 27.0 /xm and 14.0 /xm, respectively. The symbols 
— and + denote a vortex with clockwise or counterclockwise 
circulation, respectively. The black arrow with a solid line 
indicates the direction of motion of the potential, and the 
dotted arrow indicates the direction of motion of the solitons. 
The square boxes with solid lines enclose the solitons. The 
value of the phase varies from — tt (white) to tt (black). The 
amplitude and frequency of the oscillation are 3 jim. and 190/s, 
respectively, (a) shows the density and phase profiles before 
the collision. Then, a new ghost vortex pair nucleates in (b), 
coming out of the potential. The annihilation of vortices leads 
to the nucleation of solitons in (c). These solitons move away 
from the potential in (d). 
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FIG. 10: (Color online) Transfer of solitons: The density and 
phase profiles at (a) t — 20.8 ms, (b) t — 25.5 ms , and (c) t 
— 27.1 ms are shown, where the x and y dimensions of each 
image are 32.8 /xm and 34.0 /xm, respectively. The overall 
density profile at t = 59.0 ms is shown in (d), where the x 
and y dimensions of each image are 145 /xm and 34.0 /xm, 
respectively. The symbols — and + denote a vortex with 
clockwise or counterclockwise circulation, respectively. The 
value of the phase varies from — tt (white) to tt (black). The 
amplitude and frequency of the oscillation are 2 /xm and 280/s, 
respectively. The solitons, indicated by the square box with 
solid lines, nucleate by the potential in (a), moving toward the 
surface of the condensate in (b)-(c), where the dotted lines 
show the solitons. The solitons cut on the condensate in (d). 



arate and the resulting vortices move toward the bow of 
the condensate by the same mechanism in Fig .4. Note 
that these dynamics do not always occur, because the 
oscillating potential occasionally causes the solitons to 
crash, hindering the above dynamics. 

Migration of vortices - The vortices that do not change 
into solitons migrate in the condensate in Fig. 8 where 
two vortices among them are depicted by + and — sym- 
bols. First, the vortices near the surface of the conden- 
sate move towards the bow as shown in Fig. 8(a). Sec- 
ond, the vortices that have reached the bow reconnect 
with each other in Fig. 8(b). Third, the new vortex 
pair returns to the center of the condensate in Fig. 8(c). 
Figure 8 also shows that the surface is distorted by the 
phonons emitted by the oscillating potential. 

Thus, the oscillating potential leads to the synergy dy- 
namics and migration of vortices. Migration is also ob- 
served in a uniformly moving potential [10], but other 
dynamics are characteristic of an oscillating potential, 
which cannot occur for a uniformly moving potential. 
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FIG. 11: (Color online) Soliton decay in the center of the 
condensate: The density (left) and phase (right) profile at (a) 
t — 33.2 ms, (b) t — 36.3 ms, and (c) t — 42.4 ms are shown, 
where the x and y dimensions of each image are 32.8 /xm and 
14.0 /xm, respectively. The symbols — and + denote a vortex 
with clockwise or counterclockwise circulation, respectively. 
The value of the phase varies from — tt (white) to tt (black). 
The amplitude and frequency of the oscillation are 3 /xm and 
190/s, respectively, (a)-(b) show the solitons moving away 
from the potential, which are depicted by the square boxes 
with solid lines. The soliton decays to a vortex pair in (c). 



V. NUCLEATION OF SOLITONS NEAR THE 
OSCILLATING POTENTIAL 

An oscillating potential with an amplitude smaller 
than L leads to different dynamics than the case of a 
large amplitude. The most distinctive dynamics occur in 
the collision between vortex pairs and the oscillating po- 
tential. In the following, we show the dynamics related 
to the collision between the pairs and the potential, and 
the propagation and decay of solitons. 

Nucleation of solitons near the oscillating potential oc- 
curs in the case of a small amplitude, shown in Fig. 9. 
The mechanism of vortex nucleation in Fig. 9(a) is the 
same as the case for a large amplitude. The difference 
occurs in the collision between vortex pairs and the os- 
cillating potential in Fig. 9(b). The distance between 
the vortices and the potential is closer than that in the 
case for a large amplitude, which means that vortices in 
the condensate are close to the ghost vortices inside the 
potential. Therefore the annihilation of vortices occurs 
easily, leading to the nucleation of solitons in Fig. 9(c). 
Thereafter, these solitons move away from the potential 
in Fig. 9(d). 

This phenomenon is attributed to the relation between 
the nucleation time and the quarter period of the oscilla- 
tion. Generally, it takes a time ty for a moving potential 
to nucleate a vortex pair. If the time ty is larger than 
the period tp = t:/2uj^ the nucleation cannot occur since 
the velocity of the potential decreases after the period tp 
before the time ty passes. In the case of a large ampli- 
tude, tp of the oscillation is much longer than ty because 



the frequency uj of the oscillation is small. Thus, the po- 
tential can move in one direction for a long time tp > ty, 
nucleating the vortices. After the nucleation, the poten- 
tial still moves in the same direction, so that the vortices 
part from it as shown in Fig. 2(c). On the other hand, a 
potential with a small amplitude needs a large frequency 
to cause nucleation, which implies that the potential can 
move for a short time tp > ty. In this case, the vortices 
are nucleated just before the potential changes the direc- 
tion of the motion, as shown in Figs. 9(a) and (b). Thus, 
the situation in which the distance between the vortices 
and the potential is close causes annihilation of vortices 
near it, and the solitons are nucleated in Fig. 9(c). 

There are two patterns for the dynamics of solitons 
after nucleation. One is that the solitons move toward 
the surface without breaking, and the other is that they 
decay to a vortex pair in the center of the condensate. 

We show the former pattern in Fig. 10, where the 
amplitude and frequency of the oscillation are 2 /xm and 
280/s. The solitons are nucleated near the potential in 
Fig. 10(a) through the mechanism in Fig. 9. In this case, 
the solitons move away from the oscillating potential and 
toward the surface of the condensate, making cuts in it as 
shown in Fig. 10(b) and (c), where the dotted lines show 
the solitons. The oscillating potential keeps nucleating 
solitons moving toward the surface, so that several cuts 
are made in Fig. 10(d) in contrast to the case for large 
amplitude in Fig. 8(c). 

We show the latter pattern, in which the solitons decay 
to a vortex pair in the center of the condensate, in Fig. 
11, where the amplitude and frequency of the oscillation 
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FIG. 12: (Color online) Case for high velocity: The density 
profiles at (a) t = 5.89 ms, (b) t = 8.84 ms, and (c) t = 13.7 
ms are shown, where the x and y dimensions of each image 
are 28.0 /xm and 14.0 /xm, respectively. The overall density 
profile at t = 265 ms is shown in (d), where the x and y di- 
mensions of image are 145 /xm and 34.0 /xm, respectively. The 
symbols — and + denote a vortex with clockwise or counter- 
clockwise circulation, respectively. The black arrows indicate 
the direction of motion of the potential. The amplitude and 
frequency of the oscillation are 10 /xm and 90/s, respectively. 
A vortex pair is nucleated in (a) and the potential nucleates 
another pair in (b). As a result, multiple vortex pairs in (c) 
appear in the condensate. The potential continues nucleating 
multiple pairs, which leads to the distorted condensate in (d). 



9 



(a) 



e/L 



2.2 
2.0 
1.8 
1.6 
1.4 
1.2 
1.0 
0.8 
0.6 
0.4 



hig. d^b 1 ^ hig. lid 










III 



















V/Vk 0.6 



0.1 0.2 0.4 0.6 0.8 

CJO/CJOk 



(b) 

III 
■ ■ 

II 



♦ ♦ 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
CJO/CJOk 



FIG. 13: (Color online) Phase diagram of the nucleation of vortices and solitons by an oscillating potential: (a) shows the e-u) 
phase diagram of the nucleation. In (a) there are four kinds of points: no nucleation of vortices and solitons (♦), reconnection 
of vortices near the potential (t) as described in Fig. 3, nucleation of solitons near the potential (•) as described in Fig. 9 and 
nucleation of multiple vortex pairs (■) as described in Fig. 12. (b) shows the v-uj phase diagram, which is related to the e-uj 
phase diagram through v = euj. The critical velocity for the nucleation of a vortex pair and multiple vortex pair is illustrated 
by two kinds of points in (b): no nucleation of vortices and solitons (♦) and nucleation of multiple vortex pairs (■). 



are 3 jum and 190/s, respectively. As shown in Fig. 11(a) 
and (b), the solitons are nucleated near the oscillating 
potential, shown enclosed by the square boxes. Then, a 
vortex pair leaving for the bow is nucleated in the center 
of the condensate as shown in Fig. 11(c). Thereafter, 
when the pair arrives at the bow, it separates and moves 
along the surface. 



VI. NUCLEATION AND DYNAMICS OF 
MULTIPLE VORTEX PAIRS 



We show the dynamics for high velocity, which means 
that the velocity is much higher than the critical velocity 
for nucleation of vortices discussed later. The higher the 
velocity of the potential, the shorter the nucleation time. 
Thus, it is possible for the potential to make multiple 
vortex pairs in Figs. 12(a)-(c), where the amplitude and 
frequency of the oscillation are 10 /im and 90/s, respec- 
tively. After the nucleation of multiple vortex pairs, the 
potential collides with them. In this case, both reconnec- 
tion of vortices and nucleation of solitons can occur near 
the potential. As a result, the condensate is filled with 
vortices and solitons in Fig. 12(d). 

The boundary of the condensate in Fig. 12(d) shows 
the sign of the granulation of the condensate [33]. If a 
potential with high velocity would oscillate for a long 
time, the whole condensate could become granular. 



VII. PHASE DIAGRAM OF NUCLEATION OF 
VORTICES AND SOLITONS BY AN 
OSCILLATING POTENTIAL 

We now consider the nucleation of vortices and soli- 
tons by the oscillating potential and systematically in- 
vestigate the dependence of the nucleation on the pa- 
rameters. To make a phase diagram for the nucleation, 
we consider what dimensionless parameters are charac- 
teristic of the system. Our system has two features: it 
is a quantum fluid and is subject to external oscilla- 
tion. The first feature is represented by quantum cir- 
culation = h/m^ from which we define the charac- 
teristic velocity = h/mL and the characteristic fre- 
quency uo,^ = h/mL'^. The second feature is represented 
by the Strouhal number [s^ peculiar to the oscillating 
field, which is expressed by VsTs/Lg, where Fg, Tg, and 
Lg are the velocity, time, and length characteristic of 
the system. In our case, the Strouhal number is e/L for 
Vs = ecj, Ts = 1/oJ and Ls = L. Then, we create the 
phase diagram of the nucleation with three dimension- 
less parameters e/L, uj/uj^ and v/v^. Since these are 
related to each other through v = euj^ there remain two 
independent parameters. Figure 13 shows the phase di- 
agram obtained by the numerical calculation. The e-uj 
phase diagram is shown in Fig. 13(a), where there are 
four kinds of points that denote no nucleation of vortices 
and solitons (♦), nucleation of vortices and reconnection 
of vortices near the potential (t) as described in Fig. 
3, nucleation of solitons after the nucleation of vortices 
near the potential (•) as described in Fig. 9, and nu- 
cleation of multiple vortex pairs (■) as described in Fig. 
12. Note that nucleation of vortices means that the os- 
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cillating potential nucleates usual vortices, discussed in 
the next chapter in detail. The velocity v of the poten- 
tial is expressed by e a;, so that we can transform the e-uo 
phase diagram into the v-uj phase diagram in Fig. 13(b) 
to show the dependence of the critical velocity for the 
nucleation of a vortex pair and a multiple vortex pair on 

The phase diagram shows that nucleation induced by 
the oscillating potential can be roughly classified into 
three regions. In region I, vortices and solitons are never 
nucleated by the potential, which only emits phonons. In 
region II, the potential nucleates a vortex pair, which is 
described in Sees. IV and V. Besides, region II can be 
classified into two sub-regions. One is a sub-region where 
reconnection of vortex pairs occurs near the potential, as 
shown in Fig. 3. In the other sub-region solitons are 
nucleated near the potential, as shown in Fig. 9. Figure 
13(a) shows that the region for reconnection of vortices 
near the potential (▼) is wider than that for nucleation of 
solitons near the potential (•) if the amplitude e is larger 
than L. The region corresponding to e/L = 1.2 ~ 2.0 is 
mostly filled with T symbols, but we observe nucleation 
of solitons (•) for e/L = 1.6 because of the following. The 
distance between the vortices and the potential is essen- 
tial for the dynamics (•) as explained in Sec. V. When 
the amplitude and frequency of the potential is slightly 
larger than the critical values for nucleation of vortices, 
the distance is small. Therefore, nucleation of solitons 
should occur even for large amplitude if the parameters 
are close to the critical value. However, the region of this 
nucleation is narrow if the amplitude is larger than L. 
Thus, the • point with e/L = 1.6 appears. The region of 
the • symbol should survive even for a large amplitude 
(e > L) if we finely sweep the frequency and amplitude. 
In region III, multiple vortex pairs are nucleated and the 
complex dynamics of vortices and solitons as described 
in Sec. VI is induced by the potential. 

It is not the case that the dynamics induced by an os- 
cillating potential approximate those of a uniformly mov- 
ing potential if the frequency uo approaches zero. This is 
because, even though the potential oscillates with low fre- 
quency, it nucleates vortex pairs with different charges, 
and the reconnection of vortices or the nucleation of soli- 
tons near the potential not caused by the uniformly mov- 
ing potential can occur. Therefore, the dynamics under 
an oscillating potential is qualitatively different from that 
with a uniformly moving potential. 



VIII. CRITICAL VELOCITY FOR 
NUCLEATION OF VORTICES 

The mechanism of nucleation of quantized vortices is 
one of the most important issues for quantum fiuids. This 
is because the vortices, which are stable topological de- 
fects, bring remarkable changes to the system once they 
are nucleated. For example, the decay of the permanent 
fiow of superfiuid helium is considered to be strongly re- 



lated to the motion of vortices, which leads to phase slip- 
page. 

Our numerical calculations suggest that the mecha- 
nism of nucleation of vortices in a system with an oscillat- 
ing potential is different from that in a system where the 
condensate fiows uniformly without any obstacles. Nucle- 
ation in a uniform system requires excitation of phonons 
since a low density region causing divergence of the quan- 
tum pressure can be formed only by the growth of the 
amplitude of phonons, which leads to the breakdown of 
Kelvin's theorem on circulation, as explained in Sec. III. 

However, in a system with an oscillating potential, a 
different mechanism of nucleation applies. In this system, 
the potential creates a low density region inside itself, 
causing divergence of the quantum pressure. This leads 
to nucleation without any growth of phonons, as shown 
in Fig. 2 where a ghost vortex pair is nucleated inside 
the potential and moves away from it to become a usual 
vortex pair. According to the numerical calculations, an 
oscillating potential continues nucleating phonons in the 
condensate, though the phonons do not grow to create a 
low density region causing the divergence when the am- 
plitude and frequency of the oscillation is near the critical 
value for nucleation of vortices. If the parameters of the 
oscillation are much larger than the critical values, the 
phonons nucleated by the potential may form a low den- 
sity region. Therefore, our numerical calculations show 
that the critical velocity for nucleation of vortices by a 
potential is smaller than that by the growth of phonons 
in our system. 

We now consider the nucleation of a single vortex pair 
in Fig. 2, rather than multiple vortex pairs in Fig. 12, 
and define the nucleation induced by the oscillating po- 
tential. Our definition of nucleation is based on the in- 
crease of the total energy of the system and the density 
profile as shown in the following. The total energy in- 
creases by the energy of the vortices when the potential 
nucleates vortices, which is roughly estimated to be the 
order of 10~^^ J in our simulation by using the energy of 
a vortex in a uniform system [38] expressed by 

D 

e, =^no— log(1.464-), (7) 

m ^ 

where D is the size of the system, which is the Thomas- 
Fermi radius Ry in our case. Figure 14 shows the time 
evolution of the total energy Et expressed by 

Et = j ri-^'^^Wr + I V\ipfdr + ^ J IVI^dr. 

(8) 

The total energy monotonically increases with oscillation 
because phonons are excited by the oscillating potential. 
For example. Fig. 14(a) and (b) show the density profile 
with a monotonic increase of the energy. In Fig. 14 (a), 
the vortices do not sufficiently escape from the bound- 
ary layer of the width of the order ^ of the potential, 
and the increase of the total energy for the oscillation 
is smaller than 10~^^ J. After Fig. 14 (a), when the 
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FIG. 14: (Color online) The increase of the total energy in the case with e = 6 /im and uj — 73.4 /s is shown, (a)-(d) show 
the density profiles at t = 562, 580, 745, 758 ms. The symbols — and + denote a vortex with clockwise or counterclockwise 
circulation, respectively. The black arrows indicate the direction of motion of the potential. The energy graph shows that there 
is a large increase for 750-850 ms, corresponding to the energy of vortices, (a)-(b) show density profiles in which the vortices 
appear only within the boundary layer, not those distant from the potential. Vortices appear in the condensate in (c)-(d), 
which is regarded as the nucleation of vortices. 



potential changes its direction of motion, ghost vortices, 
with a direction of impulse opposite to that of the usual 
vortices are nucleated inside the potential and the usual 
vortices disappear through annihilation with the ghost 
vortex pairs. Figure 14 (b) shows the density profile just 
before the annihilation. The vortices disappear within 
the boundary layer of the potential and no solitons are 
nucleated. Then, neither reconnection in Fig. 3 nor nu- 
cleation of solitons in Fig. 9 occur. Thus, vortices appear 
and disappear within the boundary layer of the potential. 
We judge that the vortices are not yet nucleated by the 
potential in this case. On the other hand, the behavior 
of the total energy in Fig. 14 changes markedly dur- 
ing t = 750-850 ms when the vortices are nucleated by 
the oscillating potential. Figure 14 (c) shows the density 
profile for the nucleation of vortices, which shows that 
vortices sufficiently escape from the boundary layer of 
the potential. An increase of the order of 10~^^ J occurs 
twice within a period of the oscillation 2t: /uj because of 
the following. The potential moves in the x direction 
and nucleates a vortex pair during the former half period 
tt/cj. Then, it starts to move back in the —x direction 
and nucleates another pair which has an impulse opposite 
to that of the previous pair during the latter half period 
tt/cj. Thus, four vortices are nucleated in Fig. 14(d), 
which is peculiar to an oscillating potential. According 
to the above, we judge that vortices are nucleated if an 
increase of total energy of the order of 10~^^ J occurs 
twice, as shown in Fig. 14. Using this criterion for nu- 
cleation, we find that the critical velocity of the vortices 
Vc/v^ is 0.385 - 0.514 as shown in Fig. 13(b) (♦) by 
changing the amplitude and the frequency of the oscilla- 
tion by 1 Ilia and 10 /s, respectively. 



There are two candidates for the critical velocity for 
nucleation of vortices in our system. One is the sound 
velocity Cs and the other is v,^ = n/L. These veloci- 
ties are based on the Landau instability, which occurs at 
V = iam[Ee/pe]i where Ee and Pe are the energy and mo- 
mentum of an excitation, respectively. If the energy and 
momentum of a Bogoliubov excitation (a vortex pair) is 
chosen to be and respectively, we obtain Cg {v,^). 
In the following, we consider whether these velocities are 
the critical velocity in our system. 

The sound velocity Cs near the potential is often com- 
pared with the critical velocit y for nucleation of vor- 
tices in atomic BECs plol, [27l, l28l, Wt\ . However, the ve- 
locity Cs cannot be the critical velocity Vc in our system 
where a local potential moves in the condensate. The first 
reason is that the standard representation Cg = \J griQ/m 
is rewritten as Cs = k,/2\^7t^ including the coherence 
length ^ = h/ ^/2mnQg and the quantum circulation 
but the size L of the potential is not included. Sasaki 
et al. [11] have reported that the critical velocity Vc 
is dependent on L, which is obtained by the numerical 
simulation of the GP equation with a uniformly moving 
potential. The second reason is that the velocity Cg is 
just the critical velocity for Landau instability related to 
phonon emission when the condensate flows in a uniform 
system with dissipation and noise. Our simulation shows 
that the vortices can be nucleated without the growth of 
phonons, which means that the velocity Cg is not directly 
related to the nucleation of vortices, being just the critical 
velocity for which Bogoliubov phonons are excited in a 
uniform system. When the critical velocity Vc/vk in Fig. 
13(b) is renormalized by the sound velocity c^, Vc/cs is 
about 0.292 ~ 0.391. A similar result has been reported 
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by Jackson et al. The disagreement between the 

critical velocity and the sound velocity is attributable to 
this discussion. 

The velocity cannot be the critical velocity Vc either 
since the healing length ^ is not included in the velocity 
v^. The coherence length can influence the critical ve- 
locity for nucleation of vortices in atomic BECs, which is 
different from the case of superfluid helium. In helium, 
the critical velocity is considered to be almost dominated 
by the quantum circulation tv, from which the velocity is 
estimated to be [35[. This idea is based on the sit- 
uation where the size of the oscillating object is much 
larger than the coherence length in helium; the size of 
the objects is usually larger than the order of /im and 
the coherence length is of the order of A. The small 
coherence length is caused by the strong interaction be- 
tween helium atoms. However, in dilute atomic BECs, 
the coherence length is slightly smaller than the size of 
the oscillating potential because the system is dilute and 
the interaction between particles is weak. The typical co- 
herence length and size of the potential are of the order of 
0.1 /im and 1 /im, respectively. Thus, the critical velocity 
in atomic BECs does not conform with the rough esti- 
mate in which the coherence length is not taken into 
account [s^. Our numerical calculations show that the 
critical velocity Vc/v^ is 0.385-0.514 as shown in Fig. 13 
(♦), which is considerably smaller than the characteristic 
velocity = 934 /im/s. 

From the above discussions, showing the importance 
of hi^ L, and ^ for the critical velocity, we propose an ex- 
pression for the critical velocity including these variables: 



where ^ is the coherence length near the center of the 
condensate and a is a parameter. In our case, a is about 
11-19 because the critical velocity is 360-480 /im/s. This 
means that the potential is covered with a low density 
region of the width of the order of ^, and the effective 
size of the potential is larger than L. This expression ap- 
proaches and Cs for the two limiting cases. If L is much 
larger than ^, Vc becomes v^. On the other hand, Vc is 
close to the order of Cs if ^ is much larger than L, which 
means that the system becomes uniform. Note that this 
discussion should be correct for >> /i because of the 
following. If the strength Vq of the oscillating poten- 
tial satisfies Vq » /i, an increase of Vq only enlarges L 
without changing the density profile inside the potential, 
where the density is very low for >> /i. Therefore, Eq. 
(9) contains information on Vq through L. However, if 
Vq becomes comparable to /i, the condensate penetrates 
into the potential. In this case, Eq. (9) may be insuffi- 
cient since the density profile inside the potential changes 
dependently on Vq. 

In our case, the oscillating potential is treated, but 
Eq. (9) would also be valid in a system with a uniformly 
moving potential. Neely et al. ^Oj observed the critical 
velocity in a trapped system with a uniformly moving 
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FIG. 15: (Color online) Dynamical critical phenomena for 
nucleation of vortices: The relation between the nucleation 
time tv and frequency uj is shown. The amplitude of the 
oscillation e is set to be 6 /im. The time tv diverges as the 
frequency decreases to approach cjc = 73.5 /s. 



potential and found Vc = O.lcs. In their experiment, the 
healing length and size of the potential were ^ = 0.3 /im 
and L ^ 20 /im, respectively. As Eq. (9) is applied to the 
result, the parameter a is about 21. Moreover, we obtain 
a ~ 20 when Eq. (9) applies to the critical velocity ob- 
tained by Sasaki et al. [11|, which is the numerical result 
of the GP equation with a uniformly moving potential. 
Thus, a in Eq. (9) takes a similar value, a = 10-20, be- 
ing independent of whether the system has a uniformly 
moving potential or an oscillating potential. We consider 
that the value of the parameter a has the same order in 
these systems if the shape of the potential is the same 
and the strength of the potential is much larger than the 
chemical potential. 

IX. DIVERGENCE OF NUCLEATION TIME OF 
VORTICES 

We now address the issue of vortex nucleation with the 
relation of critical phenomena. Customary critical phe- 
nomena are classified into two categories, static and dy- 
namic. Static phenomena concerning equilibrium prop- 
erties mean that a system shows a characteristic power 
law for order parameters, susceptibility, specific heat, etc. 
near the critical region. Dynamic phenomena concerns 
time dependent phenomena such as relaxation time, dif- 
fusion constant, etc. Such dynamic critical phenomena 
have been reported in quantum fluid dynamics. Ther- 
mal counterflow of superfluid ^He is known to experimen- 
tally show a transition from a laminar state to turbulence 
states consisting of a tangle of quantized vortices when 
we increase the relative velocity between the normal fluid 
and superfluid [39] . When the turbulence is generated in 
a tube of circular cross section, a vortex-line density is 
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observed to develop from a low-density state (T-I) to a 
higher- density state (T-II) that can be associated with 
the homogeneous state. Griswold et al. observed that 
a large fluctuation of the vortex-line density sharply in- 
creased the time constants of the system near the transi- 
tion between T-I and T-II [40|. The transition to QT in 
thermal counterflow can be understood in the light of a 
recent numerical simulation of the vortex fllament model 
based on the full Biot-Savart law (4l|. The relaxation 
time to a turbulent state increases sharply near the crit- 
ical region between the laminar and the turbulent states 
because of the large fluctuation of the vortex-line density 

The nucleation time ty of vortices by an oscillating po- 
tential diverges near the critical amplitude and frequency 
of the oscillation. We expect that there is a critical re- 
gion where a power law between them is realized because 
Huepe et al. [28] considered a power law with saddle- 
node bifurcation in a system with a uniformly moving 
potential. The nucleation time ty is determined by the 
characteristic increase of the total energy and the appear- 
ance of vortices in the density proflle, as discussed in Fig. 

14. We show the relation between ty and the frequency 
uj with the flxed amplitude e = 6 jam in Fig. 15. We flnd 
that ty diverges as uj approaches ujc = 73.5 /s in Fig. 

15. This divergence means that the characteristic time 
becomes long when the dynamical response of the system 
drastically changes across Uc- However, we could not flnd 
a critical region in this work, probably because the region 
is too narrow to appear in our calculation. This region 
may be wide for different values of amplitude e. 
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FIG. 16: (Color online) Increase of the total energy: The 
four kinds of points denote the increase of the total energy 
for UJ = 70/s (+), UJ = 90/s (x), uj - 110/s (*), uj = 130/s 
(□) with e = 5 /im. For the + and x points, vortices are not 
nucleated, while vortices are nucleated for * and □. 



of phonons. However, in the case of nucleation of vor- 
tices {uj = 110 /s and uj = 130 /s), the energy greatly 
increases. Thus, there is a clear difference between the 
two cases, as reported by Jackson et al. [27]. 



XI. DISCUSSION 



X. HEATING OF THE CONDENSATE 

It is possible for an oscillating potential to increase the 
temperature of the system, and the GP model may break 
down. We can numerically calculate the total energy as 
shown in Fig. 16, qualitatively estimating the increase 
of the temperature by using the speciflc heat C of the 
equilibrium state for an ideal Bose gas, which is expressed 

by 

h'^UJx^y 

where ({n) is the zeta function. In the case of e = 5/im 
and u = 130 /s in Fig. 16, the increase of the total 
energy is 2.44 x 10~^^ J for 400 ms and the speciflc heat 
is 5.40 X 10~^^ J/K, which leads to a temperature increase 
of 0.45 nK. This is very small compared with the typical 
temperatures in experiments, which are well described 
by the GP equation. Therefore we conclude that the GP 
model is still valid against the heating in our simulations. 

Figure 16 also shows that the increase of total energy 
is different depending on whether the vortices are nucle- 
ated. The energy for the no vortex situation {uj = 70 
/s and cj = 90 /s) slightly increases due to nucleation 



We have performed numerical calculations with the GP 
equation under the condition that the aspect ratio of the 
shape of the condensate is flxed. One may be concerned 
about the dependence of the aspect ratio on the dynam- 
ics. However, we believe that the dependence is weak. 
The dynamics induced by an oscillating potential can be 
classifled into two kinds. One is the dynamics near the 
potential, which includes the reconnection of vortices and 
nucleation of solitons near the potential. These dynam- 
ics are barely concerned with the surface, occurring with 
an arbitrary aspect ratio. The other is the dynamics of 
vortices and solitons related to the surface of the conden- 
sate, including the migration of vortices in Fig. 8 and the 
collision and collapse of solitons in Fig. 7. The migra- 
tion is considered to have a weak dependence because 
the vortices only move along the surface. On the other 
hand, the collision and collapse can have a dependence 
on the aspect ratio. The direction of motion of the soli- 
tons depends on the curvature of the surface, so that the 
collision and collapse may not occur. 

Our numerical calculations do not contain noise in the 
initial state, so that all our results have mirror symme- 
try about the x axis. However, our results should be 
consistent with experiments with noise and fluctuation 
breaking the symmetry of the system. Our results for 
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the dynamics of vortices and solitons are dominated by 
their behavior near the potential, such as the reconnec- 
tion of vortices and nucleation of sohtons near the po- 
tential. These dynamics are local, so that the noise do 
not influence them. Though the noise would break the 
mirror symmetry of the condensate, the overall dynamics 
would not be qualitatively affected. 

The results in this paper are based on the two- 
dimensional (2D) GP equation, from which we can in- 
fer the dynamics of vortices and solitons for the three- 
dimensional (3D) case. A vortex pair and soliton with 
linear configuration in 2D correspond to a vortex ring 
and a soliton with planar configuration in the 3D case. 
Thus, the dynamics in 3D would be obtained by replac- 
ing the vortex pair and linear soliton with a vortex ring 
and planar soliton . An oscillating potential in 3D can 
create a vortex tangle in the condensate and lead to QT. 

Finally, we consider the appropriate parameters for an 
oscillating potential for the growth of QT. QT needs vor- 
tices and hence the regions II and III in Fig. 13 would 
be appropriate. However, region II with nucleation of 
solitons near the potential will not develop to QT as the 
number of the nucleated vortices is insufficient. There- 
fore, we suggest region II with a large amplitude com- 
pared with the size of the potential and region III as 
appropriate parameters for the growth of QT. An oscil- 
lating potential with the parameters in regions II and III 
may lead to granulation of the condensate. In fact, the 
sign appears in Fig. 12 (d), where the boundary of the 
condensate is granular. 

XII. CONCLUSION 

We numerically calculated the two dimensional GP 
equation to investigate the dependence of the parame- 
ters on the dynamics of the condensate induced by an 
oscillating potential. We roughly obtained four kinds of 
dynamics. First, the oscillating potential does not nucle- 
ate vortices and solitons; only the emission of phonons 



occurs. Second, the potential nucleates vortices and re- 
connection of vortices near the potential in Fig. 3 oc- 
curs. Thereafter, the vortices migrate in the conden- 
sate or transform into the solitons as shown in Sec. IV. 
Thirdly, the potential nucleates solitons near itself in Fig. 
9, which is very different from the case in Sec. IV. Af- 
ter nucleation, the solitons transform into vortex pairs or 
transfer to the surface of the condensate as shown in Sec. 
V. Fourthly, multiple vortex pairs are nucleated by the 
potential, which forms the distorted condensate in Fig. 
12. These dynamics depend on the parameters of the 
oscillating potential, which is represented by the phase 
diagram in Fig. 13. 

We minutely discuss the nucleation of vortices induced 
by an oscillating potential. The key to the nucleation is 
the divergence of the quantum pressure, which leads to 
the breakdown of Kelvin's theorem on circulation. Thus, 
nucleation requires a low density region to cause the di- 
vergence, from which it follows that the mechanism of the 
nucleation has two processes. One is the process in which 
the vortices are nucleated through the growth of the am- 
plitude of phonons. The other is that the low density 
region inside the potential causing the divergence leads 
to nucleation of vortices without the growth of phonons. 
We propose a new expression of Eq. (9) for the critical 
velocity Vc of the nucleation, including both the coher- 
ence length ^ and the size L of the potential. 

We treated the nucleation of vortices by an oscillat- 
ing potential as a problem of dynamical critical phenom- 
ena, and investigated the relation between the nucleation 
time and the frequency. Our numerical calculation re- 
vealed that the nucleation time diverges near the critical 
frequency, but we unfortunately could not confirm the 
power law between the time and the frequency. 
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